#include "TH1D.h"
#include <vector>
#include <iostream>
#include "TROOT.h"
#include "TMath.h"
#include "TCanvas.h"
#include "TRandom.h"
#include "TF1.h"
#include "TFile.h"
#include "TChain.h"
#include "TClonesArray.h"
#include "TLorentzVector.h"
#include <TStopwatch.h>
#include "TString.h"

using namespace std;
// ******** varialbe for mathmatics ***************
const double pi = TMath::Pi();
const double PI = TMath::Pi();




// ******** Flow setting Parameters ***************

const int NH = 4;
const double input_vn[NH] = {100, 0, 0.06, 0.01};


enum NdetType{kTargetDetector, kEP4Det, kEP8Det, kEP16Det, k32Det, k64Det,kIealDet, kNDetector};
const int ND = kNDetector;
const double det_resol[ND]={0, pi/2, pi/4, pi/8, pi/16, pi/32, 0};
TString Title[ND]={"kTargetDetector" , "kEP4Det", "kEP8Det", "kEP16Det", "k32Det", "k64Det","kIdealDet"};

